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Abstract 

We draw some rigorous conclusions about the functional properties of the 
yU — p relation in the Hubbard model based on symmetry considerations 
and unitary transformations. It is shown that the charge susceptibility 
reaches its local extreme at half-filling. Exact expressions are obtained in 
two limiting cases. 
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1 Introduction 



Exact results are of great interests in discussions of strong correlations [T[]. The important 
one-dimensional strongly correlated models such as Hubbard model and Kondo model can be 
exactly solved by means of Bethe ansatzp[||^. Nevertheless, exact knowledge of higher dimen- 
sional models are very rare. Remarkably, the 77-paring mechanism^ can be employed to discuss 
the possible existence of the off-diagonal long range order(ODLRO) at both zero and finite 
temperatures |0 . The exact knowledge about the chemical-potential relation and magnetiza- 
tion , which are important aspects of the statistics, is still called for. Because of this situation, a 
number of approximate approaches are applied, as well as computational simulations. It should 
be realized that every reliable approximation must coincide with the exact one qualitatively. 
That is, it should have the properties of the exact solution. Therefore, the exact properties of 
the /i — p relation should be studied. This is the motivation of this paper. The layout of this 
paper is as follows. Section 2 discusses the general exact properties of /i — p relation. The last 
section is devoted to conclusional remarkes. 



2 Properties of the fi — p relation 

The Hamiltonian of the generic Hubbard model on lattice A is 

H = Y.Y1 tijctcja + UY^ rii^n.i (1) 

ij cr i 

and the grand canonical partition function is 

Z{ij,(3;t,U) = TTexp[-(3{H - 12N)] (2) 

where and Qo- are the creation and annihilation operators of the electrons with spin a =|, | 
at site i. The hopping matrix {tij} is required to be real and symmetric. The number operators 
are njo- = cj^Qo-, while the U denotes the on-site Coulomb interaction. It is further assumed 
that the the lattice A is bipartite in the sense that it can be divided into sublattices A and 
B such that tij = whenever {ij} G A or {ij} E 'B. N = J^i.o-'^io- is the number operator of 
the electrons. Though the desired fi — p relation is of the form /i = fi{p, P) in statistics, the 
properties to be discussed below is that of the inverse relation p = p{p,P). 





We first derive a recursion relation concerning the moments of tlie number operators. For 
tfiis purpose, we make use of tlie complete particle-hole transformation 

PC,,P-' = e{z)cl (3) 

where e{i) = 1/ — lifie A/B. Under this transformation, the grand canonical Hamiltonian 
undergoes the following transformation 

P(H - iiN)p-^ + U(Na -N)- ii(2Na - N) (4) 

where A^a is the total number of sites of the lattice A. Therefore, the partition function satisfies 
the following relation 

Z{fi, 13- 1, U) = exp[-/3(t/ - 2fi)NA]Z{U - fi, /3; t, U) (5) 

We have accordingly for n = 1, 2, 3, • • • 

< AT" >^ = Z-^ Tr PN"P-^ ex.p[-(3P(Ho - iiN)p-'^] (6) 
= Z-^Tr(2iVA-iV)"exp[-/5P(i/o-/^iV)p-^] (7) 
= <{2NA-Nr>u-i, (8) 

This relation holds also locally. Consider the local density operator = Z^o-cl^Qo-, using the 
transformation (3), the same reasoning leads to 

< >^=< (2 - >u-^ (9) 

This local relation can not be obtained from (8) directly. In particular, we have at half-filling 
for which /j, — — U/2 from (9) < rii >^^— 1. The global relation (8) gives then 

< TV" >^,=< (27Va - TV)" (10) 

That is 

n 

[1 + (-1)"+^] < TV" >^,= ^Ci(2iVA)'(-l)"-' < iV"-^ (11) 

1=1 

where CL is the combinatorial numberiCf, = The first three cases for odd n are 

<N>^^ = TVa (12) 
<7V3>^^ = 3iVA < iV2 -2iVl (13) 
<7V^>^, = 5Na<N''>^^-20NI<N^>^^+16NI (14) 



In general, eq.(lO) turn out to be identities for even n and it provides non-trivial information 
only for odd n. On the other hand, since the fluctuation 

is of order < iV >, we may introduce a function a via 



^— < ^ >M=< N^>^.-<N>1 (15) 



^§-^<N >/.= <N>^ (16) 

One can see that the derivatives of < > with respect to /x exist at every order. Obviously, 
one has 

a > (17) 
The function a is related to the compressibility which is defined by 

1 dV 

- = -V^dp)T,N (18) 

as follows. Since p = N/V, we may express k, as 

I, dp 1 dp dp 

Using the thermodynamic relations 

/^=(^)P,T, V=i—)r,N (20) 



where $ is the thermodynamic potential, we have 

.dp. .dV. 



(21) 



therefore 



1 dp a/3 

«=— (t^)t = — (22) 
op p 

This relation is more thermodynamical than statistical. It holds for both fermi liquids and 
non-fermi liquids. For the electron gas at zero temperature, this relation can be derived in 
another way0. In the Hubbard model, the density is defined by p =< N > /N^. Since < N > 
has maximum value 2N\ and is a monotonically increasing function of p, we have 

lim a = (23) 



A 



We now prove some more properties of a. Firstly, for a homogeneous system, we have 
< AT Ni^<ni> and < A^^ >= Eij < niUj >= A^a < niUQ >. Thus (15) imphes that 

1^ = 4^^[l^ 12 < ^i^o > - <no >^] (24) 

which provides another way to evaluate the compressibility. Define the density-density corre- 
lation function 

CiaJa'ifJ') =< ^i^^ja' > - < >< n-^> > (25) 

the relation 

c,.,jA^^) = c,^,Au - 1^) (26) 

can be obtained (not only for a homogeneous system) by the local relation(9). Secondly, 
applying the complete particle- hole transformation to eq.(9), we have 

< {2Nk - Nf >u-i, -<2Na-N >l_= a{n) < 2Na - N >u-^ (27) 

Expanding this relation and using eq(15) for //—>[/ — one can obtain 

[«(//) + a{U -fJ,)]<N >u-^^^ a{p)2NK (28) 

Thirdly, differentiate (28) with respect to we have 

d 

[a\ij) + a\v){-l)\,=u-^ < N >u-^, +[a{fj,) + a{U - /x)]— < N (-1)|,=[,_^ = a;'(/x)27VA 

(29) 

Therefore, at /i — /ic 

2a(/ie)(-l)^ < N >\^ = a'(/Xe)2iVA (30) 
Hence we obtain another relation 

Oi'i/J'c) = -PcY^il^c) (31) 

Using this relation, one can prove that the charge susceptibility Xc '■= ^ takes its extreme at 
half-filling, i,e, 

X'M = (32) 

This is a rigorous result and is valid in all cases at finite temperature. Unfortunately, we can 
not know whether it is a minimum or a maximum at the present time. If it is a minimum and 



the value XcifJ'c) = 0, the system will demonstrate a Mott transition. 

We next study the Taylor expansion of p = around /i — /ic- From eq(8) we have 

p(/.)=2-p(t/-/x) (33) 

If we define 

/(/i)=p(/xe + /x)-l (34) 

then 

/(//) = -fi-i,) (35) 

i.e. /(/x) is an odd function of /i. For finite temperature, p = p{fi) is an analytic function and 
can be expanded in powers oi /i — /ic, thereby, we may infer that the Taylor expansion of p(/x) 
must be of the form at finite temperature 

oo J2k+1)( \ 

+ (36) 

The average values < A'"'^ for both odd and even n must be known in order to know the 
Taylor coefficients in (36). Unfortunately, the recursion relation eq(8) can not provide us with 
useful information about < for even n . Otherwise, the exact form of — p relation 

can be immediately obtained after the sum is implemented. In any case, equation of motion is 
needed in order to calculate < N'' > for the moment. 

Though the general exact p — p relation is not accessible so far, the specific expressions for 
two limiting cases: non-interacting case and atomic case are obtainable nevertheless. In the 
first case, U — 0, if the lattice is periodic, we can work in the k-space and we have 

where i(k) is the Fourier transform of tij. Therefore the /J, — p relation is 

P = ^ E < > (38) 
In the second case, t^j = 0, the partition function can be calculated 

Z{(3, At; C/) = Tr exp[-/?(t/ ^ ^n^a - A*^)] = [1 + 2e^'^ + e'^^^'^-^)]^^ (39) 

i 

Therefore, the average density is 



This relation is also exact for the case of infinite hopping range [pHI] . It can be seen that the 
Taylor coefficients must be in general both tij and U dependent as can be seen from the two 
special cases (38) and (40). The inverse relation at finite temperature = /i(p, /5) can be 
obtained. For p 7^ 2 we have 

p = i ln[p - 1 + ^(p-l)2 + p(2-p)e~2/3Me] _ 1 in(2 - p) + 2p, (41) 

while for p = 2, we can obtain only p = 00 at finite temperature otherwise the denorminator in 
(40) is finite and can be multiplied by both sides and in this way we will reach the unacceptable 
conclusion e^'^ = — 1. It should be noted that discontinuities may appear at zero temperature. 
For instance, if pc > we have from (40) 

[ hm^^o, i ln[|Me-2/^^i + 2p, = if p<l 
jimp=K, z/p=l 
[ 2pc if P> ^ 

while for Pc < and p 7^ 2, lim^j^oo P = Pc- 

For tij 7^ 0, numerical and analytical results show that there is a critical value Uc such that 
for U > f/c at T = 0, there is a gap in the density of states and a discontinuity of p[0]. 



Introducing the dimensionless parametres t = pt,fi = Pp,U = PU, p is a function of these 
parametres: p = p{t,p, U). Since under the Lieb-Wu transform in 0, t ^ — t, the dependence 
of t is in fact through . Hence we have 

piP, p; t, U) = p(p; t» = 1 + E (o^l^^i ' '^^) ^^2) 

fc=0 l^/C+ij! 

3 Discussions 

We discussed exact properties of the chemical potential-density relation by means of the 
particle-hole transformations and the symmetries of the Hubbard model. Similar discussion 
applies to Kondo lattice model and periodic Anderson models as in These exact features 
can help us to evaluate the reliabilities of the various conclusions from different approxima- 
tions. Yet, these properties are still not enough to determine the relations completely. It can 
be inferred from the exact solutions in the two limiting cases that the general exact answer 
must be extremely complicated. 
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